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Abstract 


We study natural supersymmetry in the Generalized Minimal Supergravity (GmSUGRA). 
For the parameter space with low energy electroweak fine-tuning measures less than 50, we are 
left with only the Z-pole, Higgs-pole and Higgsino LSP scenarios for dark matter (DM). We 
perform the focused scans for such parameter space and find that it satisfies various phenomeno- 
logical constraints and is compatible with the current direct detection bound on neutralino DM 
reported by the LUX experiment. Such parameter space also has solutions with correct DM 
relic density besides the solutions with DM relic density smaller or larger than 50 WMAP9 
bounds. We present five benchmark points as examples. In these benchmark points, gluino 
and the first two generations of squarks are heavier than 2 TeV, stop ti» are in the mass range 
[1,2] TeV, while sleptons are lighter than 1 TeV. Some part of the parameter space can explain 
the muon anomalous magnetic moment within 3g as well. We also perform the collider study 
of such solutions by implementing and comparing with relevant studies done by the ATLAS 
and CMS Collaborations. We find that the points with Higgsino dominant ¥$/X} mass up 
to 300 GeV are excluded in Z-pole scenario while for Higgs-pole scenario, the points with ¥8 
mass up to 460 GeV are excluded. We also notice that the Higgsino LSP points in our present 
scans are beyond the reach of present LHC searches. Next, we show that for both the Z-pole 
and Higgs-pole scenarios, the points with electroweak fine-tuning measure around 20 do still 
survive. 
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1 Introduction 


Undoubtedly, the gauge coupling unification of the strong, weak and electromagnetic interac- 
tions of the fundamental particles is a great triumph of the supersymmetric (SUSY) version 
of the Standard Model (SM) of particle physics [1], which henceforth will be called as Super- 
symetric SM (SSM). The SSM predicts the existence of SUSY partners of all the known SM 
particles. Interestingly, the existance of these particles can help us to understand the stabiliza- 
tion of the electroweak (EW) scale and thus solves yet another daunting problem of particle 
physics named as the gauge hierarchy problem [2]. In addition, the Minimal SSM (MSSM) also 
predicts the Higgs boson mass (mp) should be smaller than 135 GeV [3]. Indeed, the ATLAS 
and CMS Collaborations of the Large Hadron Collider (LHC) have discovered a SM-like Higgs 
boson h with mass mp= 125 GeV [4, 5]. This adds yet another feather in the hat of the SSM. 
The SSM also predicts that with R-parity conservation, the Lightest Supersymmetric Particle 
(LSP) such as neutralino is an excellent dark matter candidate [6, 7]. And the electroweak 
symmetry can be broken radiatively due to large top quark Yukawa coupling, etc. All these 
observations give us some hints that we are on the right track. 

The existence of the SM-like Higgs boson with mass mp, ~ 125 GeV requires the multi- 
TeV top squarks with small mixing or TeV-scale top squarks with large mixing. This raises a 
question on the naturalness of the MSSM and generates the fine-tuning problem. However, the 
null results of the LHC-Run2 and the ongoing LHC SUSY-searches have not found any SUSY 
evidences yet. In recent studies, the bounds on squark masses mg 2 1600 GeV [8] and gluino 
mass mg 2 2000 GeV [8] have been reported by the ATLAS and CMS Collaborations at the 
13 TeV LHC with 36 fb~! of data. This situation has put the promises of the MSSM under 
pressure. It is interesting to note that despite the SM-like Higgs mass being relatively heavy, 
there are some studies [9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21] which suggest that the 
naturalness problem in the MSSM can be solved successfully. In particular, in an interesting 
scenario, which is called as Super-Natural SUSY [17, 22], it can be shown that no residual 
electroweak fine-tuning (EWFT) left in the MSSM if we employ the No-Scale supergravity 
boundary conditions [23] and Giudice-Masiero (GM) mechanism [24] despite having relatively 
heavy spectra. Some people might think that the Super-Natural SUSY might have a problem 
related to the higgsino mass parameter u, which is generated by the GM mechanism and is 
proportional to the universal gaugino mass Mj/2, since the ratio Mj/2/ is of order one but 
cannot be determined as an exact number. This problem, if it is, can be addressed in the 
M-theory inspired the Next to MSSM (NMSSM) [25]. Also, see [26], for more recent works 
related to naturalness within and beyond the MSSM. 


In order to quantify the amount of fine-tuning (FT), we need to define the fine-tuning 
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measures. In literatures, we can find the high energy fine-tuning measure Aggyz—paq defined 
by Ellis, Enqvist, Nanopoulos and Zwirner [27], as well as Barbieri and Giudice [28], and the 
high energy and electroweak fine-tuning measures Apg and Agw defined by Baer, Barger, 
Huang, Michelson, Mustafayev and Tata [29, 30]. Usually, we have Agw < Ago < Ang. One 
can show that Agw ~ Agg for some scenarios [31]. 

This work is a continuation of our phenomenological studies of Generalized Minimal Su- 
pergravity Model (GmSUGRA) [32]. In Refs. [33, 34], we showed that in GmSUGRA, we 
have varieties of dark matter scenarios such as A-resonance, Higgs-resonance, Z-resonance, 
stau-neutralino coannihilation, tau sneutrino-neutralino coannihilation compatible with vari- 
ous phenomenological constraints. In addition, we showed that the Higgs coupling and muon 
anomalous magnetic moment measurements can constrain the parameter space effectively. In 
this work, we concentrate on the dark matter solutions which not only have low EWFT (that 
is Agw < 50), but also are consistent with current direct detection bounds reported by the 
LUX Collaboration [35]. In our scans, we find that the light stau-neutralino coannihilation 
points do not satisfy Agw < 50. Also, the Higgsino LSP points are still natural and viable, but 
they cannot be probed at the current LHC searches. We find that only Higgs-pole and Z-pole 
solutions fulfil the above mentioned criteria. Therefore, we will only consider these two type 
of resonance points in more details. In these two scenarios, a subset of solutions satisfy the 5c 
dark matter relic density WMAP9 bounds while the other solutions have relic density beyond 
the 50 bounds. We present five benchmark points as examples of the parameter space under 
consideration, where one of them has the Higgsino LSP. In these benchmark points, gluino and 
the first two generations of squarks are heavier than 2 TeV, top squars ti» are in the mass 
range [1,2] GeV, while sleptons are less than 1 TeV. Some part of the parameter space can 
also explain the muon anomalous magnetic moment within 30 [36]. Furthermore, we consider 
the constraints on such solutions from the direct searches for the SUSY particles at the LHC. 
In order to realize small fine-tuning and satisfy experimental constraints simultaneously, only 
electroweakinos (neutralinos and charginos) and stau are light and could be explored at the 
current LHC searches. We study various electroweak Drell-Yan production processes where one 
could produce neutralinos which could decay through on-shell or off-shell Z™ (x? > Z™¥®) 
or h*) (X? + hZ). We will give more details about the our analyses later in this paper. 
We display various plots showing that the relevance of different decay modes depends on mass 
spectra and will significantly influence collider searches for these particles. The dominant decay 
channel of ¥§ for samples of Z-pole is ¥8 > LZ when the mass difference meg —mM<o is small. 
Once the decay into Higgs boson is kinematically possible, branching ratio to ¥9h increase with 
increasing of mgg — m 


zo and become the dominant channel when myo — mgo Z 140 GeV. The 
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decay channels of %¥ is always Y¥f — eW. We also find that for our present work, 31+ Py 
and 21+ Ér give the best sensitivity at the LHC searches where electroweakinos decays to 
multi-leptons. We use suitable kinematic variables to discriminate signals from backgrounds. 
We show the 95% C.L. exclusion results of the LHC electroweakinos searches in the Mxo-Mso 
plane and m o-Apw plane. It can be seen from these plots that higgsino dominant X3/XP with 


mass up to 300 GeV are excluded in case of Z-pole while for Higgs-pole scenario, points with 
Xe mass up to 460 GeV are excluded. Moreover, it can also be noticed that Z-pole solutions 
with small Agw are easy to be explored, whereas solutions with large Agw are hard to exclude 
but for the Higgs-pole, many points with Agw up to 50 could by excluded by electroweakino 
searches with tau final states. Finally, we notice that for both the Z-pole and the Higgs-pole, 
samples with Agw ~ 20 could still survive, indicating naturalness of this SUSY framework. 
The remainder of this paper is organized as follows. We present our model in Section 2. 
We discuss EWFT measure in Section 3. Section 4 is devoted for scanning procedure and 
phenomenological constraints. Our results for focused scans are shown in Section 5 while 
results for the LHC searches are presented in Section 6. A summary and conclusion are given 


in Section 7. 


2 The Electroweak SUSY from the GmSUGRA in the 
MSSM 


In GmSUGRA, at the GUT-scale, we can write the generalized gauge coupling relation and the 


generalized gaugino mass relation as follows 


1 1 1 1 
~-i =K(A--), (1) 
Q2 Q3 Qı Q3 
M.: M. M. M; 
BB = (S-S) (2) 
Q2 Q3 Qı Q3 


where k is the index of these relations since it is invariant under one-loop Renormalization 
Group Equation (RGE) running. For more details about the model, please see [32]. 

Another important feature of GmSUGRA is that we can realize Electroweak SUSY (EW- 
SUSY). In this scenario, we can have the sleptons and electroweakinos within one TeV while 
squarks and/or gluinos can be in several TeV mass ranges [37]. Assuming gauge coupling uni- 
fication at the GUT scale (a; = a2 = a3) and using k = 5/3, we obtain a simple gaugino mass 
relation from Eq. (2) 

My — M; = 3 (Mi — Ms) (3) 
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It is straightforward to notice that the universal gaugino mass relation Mı = My = M in 
the mSUGRA, is just a special case of this general one. This is why we call it Generalized 
mSUGRA. We will choose M, and Mp, to be free input parameters, which vary around several 
hundred GeV for the EWSUSY. We can now write Eq. (3) for M3 as: 


5 3 
2 — — — 4 
M3 2 Mı 2 Mə , ( ) 


which could be as large as several TeV or as small as several hundred GeV, depending on 
specific values of Mı and Mp. 

The general SUSY breaking (SSB) soft scalar masses at the GUT scale are given in Ref. [38]. 
Taking the slepton masses as free parameters, we obtain the following squark masses in the 
SU(5) model with an adjoint Higgs field 


5 1 

mi, = mg) + im, (5) 
5 2 

Mize = 3 (mo)? — SMe > (6) 
5 2 

mbe _ zmo)? — 3M, > (7) 


where mg, Mec, Mpe, Mz, and Mge represent the scalar masses of the left-handed squark 
doublets, right-handed up-type squarks, right-handed down-type squarks, left-handed slep- 
tons, and right-handed sleptons, respectively, while m¥ is the universal scalar mass, as in the 
mSUGRA. In the EWSUSY, m; and mge are both within 1 TeV, resulting in light sleptons. 
Especially, in the limit my > mz JĀ WE have the approximated relations for squark masses: 

2 2 


2m% ~ mM; ~ Mhe In addition, the Higgs soft masses mg, and m,,, and the trilinear soft 


terms Ay, Ap and Ag can all be free parameters from the GmSUGRA [37, 38]. 


3 The Electroweak Fine Tuning 


As we mentioned earlier that in this work we are interested in solutions with low EWFT. We use 
the (7.85) version of ISAJET [39] to calculate the FT conditions at the EW scale Mew. After 
including the one-loop effective potential contributions to the tree-level MSSM Higgs potential, 
the Z-boson mass Mz is given by 

M3 (mint Ea) — (my, + Uy)tan’ Bo, 


ZA a ee 8 
2 tan? 8 — 1 n (8) 


where ©“ and X£ are the contributions coming from the one-loop effective potential defined 
in Ref. [30] and tan = =. All parameters in Eq. (8) are defined at the Mgw. In order to 


4 
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measure the EWFT condition we follow [30] and use the following definitions 
Cu, = [Min / (tan? 8 —1)|, Cu, = | — mz, tan? B/(tan? 8 — 1)|, Ca = |- wv’, (9) 


with each Cruden) less than some characteristic value of order M2. Here, r labels the SM and 
SUSY particles that contribute to the one-loop Higgs potential. For the fine-tuning measure 
we define 

Arw = max(C,)/(M3/2) . (10) 


Note that Agw only depends on the weak-scale parameters of the SSMs, and then is fixed by 
the particle spectra. Hence, it is independent of how the SUSY particle masses arise. Lower 
values of Agw corresponds to less fine tuning, for example, Agw = 50 implies Arw = 2% fine 
tuning. In addition to Agw, ISAJET also calculates Aps which is a measure of fine-tuning at 
the High Scale (HS) like the GUT scale in our model [30]. The HS fine-tuning measure A pg is 
given as follows 

Aug = max(B;)/(M3/2) . (11) 


For definition of B; and more details, please see Ref. [30]. 


4 Scanning Procedure and Phenomenological Constraints 


We employ the ISAJET 7.85 package [39] to perform the focused scans using parameters given 
in Section 2 to explore the parameter space having Z-resonance and Higgs-resonance solutions. 
In this work, we will focus on the solutions with relatively small EWFT Agw < 50. For full 
ranges of the parameter see [33]. 

In ISAJET, the weak scale values of the gauge and third generation Yukawa couplings are 
evolved to Maur via the MSSM renormalization group equations (RGEs) in the DR regulariza- 
tion scheme. We do not strictly enforce the unification condition g3 = gı = go at Mayr, since 
a few percent deviation from unification can be assigned to the unknown GUT-scale threshold 
corrections [40]. With the boundary conditions given at Maur, all the SSB parameters, along 
with the gauge and Yukawa couplings, are evolved back to the weak scale Mz. 

In evaluating Yukawa couplings, the SUSY threshold corrections [41] are taken into account 
at the common scale Msusy = „Mi, Mir The entire parameter set is iteratively run between 
Mz and Maur using the full two-loop RGEs until a stable solution is obtained. To better 
account for the leading-log corrections, one-loop step-beta functions are adopted for gauge and 
Yukawa couplings, and the SSB parameters m; are extracted from RGEs at appropriate scales 


mM; = m;(m;). The RGE-improved one-loop effective potential is minimized at an optimized 
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scale Mgysy, which effectively accounts for the leading two-loop corrections. The full one-loop 


radiative corrections are incorporated for all sparticles. 


It should be noted that the requirement of radiative electroweak symmetry breaking (REWSB) [42] 


puts an important theoretical constraint on parameter space. Another important constraint 
comes from limits on the cosmological abundance of stable charged particle [43]. This excludes 
regions in the parameter space where charged SUSY particles, such as 7, or t;, become the 
LSP. We accept only those solutions for which one of the neutralinos is the LSP. 

Also, we consider u > 0 and use m; = 173.3 GeV [44]. Note that our results are not too 
sensitive to one or two sigma variations in the value of m [45]. We use mP®(Mz) = 2.83 
GeV as well which is hard-coded into ISAJET. Also, we will use the notations A;, A,, A, for 
Ay, Ap and Ag, receptively. 

In scanning the parameter space, we employ the Metropolis-Hastings algorithm as described 
in [46]. The data points collected all satisfy the requirement of REWSB, with the neutralino 
being the LSP. After collecting the data, we require the following bounds (inspired by the LEP2 
experiment) on sparticle masses. 

(1) LEP2 constraints 


We employ the LEP2 bounds on sparticle masses 


Me, MG = 100 GeV, 
mz > 105GeV, 
mys > 103GeV. (12) 


(2) Higgs mass constraints 
The combined value of Higgs mass reported by the ATLAS and CMS Collaborations is [47] 


Mp = 125.09 + 0.21(stat.) + 0.11(syst.) GeV. (13) 


Due to the theoretical uncertainty in the Higgs mass calculations in the MSSM [48], we use the 


following Higgs mass bound 
122 GeV < mp, < 128GeV. (14) 


(3) LHC constraints 
We demand [8] 


mg > 2000GeV , (15) 
mg > 2000 GeV . 


(4) B-physics constraints 
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We use the IsaTools package [49, 50] and implement the following B-physics constraints 


1.6 x 107° < BR(B, > pty) < 4.2 x 107° (20) [51], (16) 
2.99 x 1074 < BR(b > sy) < 3.87 x 10°* (20) [52], (17) 
0.70 x 1074 < BR(B, > Tv) < 1.5 x 1074 (20) [52]. (18) 


(5) Electroweak Fine-Tuning constraint 
Because we consider the natural SUSY, the following constraint on fine-tuning measure 


Apw is applied 
Arw < 50. (19) 


(6) WMAP constraint 
We apply the WMAP9 bounds with 5c variation on DM density [53] 


0.0913 < Oh? < 0.1363 . (20) 


5 Results of focused scans 


We present results of focused scans in Fig. 1. In the top right and left panels we display plots in 
Apgw Vs. u respectively for Z-pole and for Higgs-pole scenarios while rescaled spin-independent 
(€0°!(x,p)) rate vs. LSP neutralino mass mgo is shown in bottom panel. Aqua points satisfy 
the REWSB and LSP neutralino conditions. Red, blue and green points represent the sets of 
points respectively with DM relic density consistent with, greater than, and smaller than 50 
WMAP9 bounds, as well as consistent with upper bounds reported by the LUX experiment. 
These points all also satisfy the bounds given in Section 4. We see that green points both 
for Z-pole and Higgs-pole scenarios have Agw in the range 20 to 50, while u is in the range 
of [100,450] GeV. Blue and red points have Agw as small as 24 and p is confined between 
[300, 450] GeV. In the bottom left panel, we show a plot in the same plane for Higgsino LSP 
solutions. Because such solutions have small relic densities, all the points are green. One can 
see that such solutions have Agw values from 20 to 50 with corresponding u values varying 
roughly from 250 GeV to 450 GeV. We will talk about Higgsino LSP solutions more with 
reference to the plots in o° (x, p)-mzo plane. 

In £0 (y, p)-mxo plot, solid black and red lines respectively represent the current LUX [35] 
and XENONIT [54] bounds. The dashed green and red lines display projection of XENON1T 
for next two years and XENONnT (total exposure of 20 t.y) [55], respectively. The factor 


€ = Oh?/0.12 for green points which accounts for a possible depleted local abundance of 
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Figure 1: Agyw vs. js for the Z-pole (left), for Higgs-pole (right) and Higgsino LSP (bottom left) 
scenario. Rescaled spin-independent (€0°!(x,p)) rate vs. LSP neutralino mass mgo (bottom 
right). Aqua points satisfy the REWSB and LSP neutralino conditions. Red, blue and green 
solutions represent the sets of points with relic density consistent with, greater than and smaller 
than 50 WMAP9 bounds, respectively. These points also satisfy the bounds indicated in 
Section 4. 
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Mio 
m 


dy, 
Mio 
Miz 
MF 2 

os1(pb) 


osp(pb) 
Qepmh? 


pkk 


BMP 1 


326 
29 
1691 
4.17 x 1071? 


45.9, 326 
337,712 
333, 704 


2374, 2542 
1173, 1731 
2375, 2561 
1717, 2433 
978 
935 
984, 909 
816, 941 
8.05 x 107+! 
1.19 x 1075 
0.106 


MMP 2 
1424 
1323.3 
1770.3 
1701.7 
851.8 
607 
98.19 
751.8 
-882.22 
3689 
-413.5 
18.9 
2562 
2697 
254 
24 
1597 
5.59 x 10710 


45, 255 
266, 658 
260, 651 
2025 
2216, 2411 
1000, 1542 
2218, 2434 
1525, 2287 
878 
821 
883, 839 
719, 828 
1.64 x 10 7 
3.38 x 10-5 
0.017 


BMP 3 
1537 
1404.1 
1981.3 
1945.9 
475.7 
132.2 
132.6 
1023 
-1203 
4981 
-238.2 
21.3 
3231 
3284 
351 
35 
2552 
6.34 x 10-10 


62, 355 
363, 882 
362, 876 


2676 
2752, 2975 
1069, 1811 
2753, 3016 
1812, 2777 

670 

032 

683 

264, 549 
7.33 x 107" 
9.07 x 1076 

0.103 


BMP 4 
1367 
1249.4 
1760.5 
1717.4 
497.6 
151.3 
131.1 
1105 
-1329.8 
5076 
-186.9 
19.8 
3306 
3203 
276 
31 
2660 
4.93 x 1071? 


62, 283 
287, 953 
286, 946 

2918 

2873, 3026 

1046, 1960 

2875, 3047 

1969, 2831 

774 
679 
786, 522 
162, 693 
1.64 x 10” 
2.40 x 10-5 
0.022 


BMP 5 
1053 
1021.8 
1169.3 
1189.1 
806.8 
849.1 
857.7 
706.8 
1084.1 
857.7 
-2915 
14.5 
2558 
487.6 
262 
29 
1597 
3.61 x 107! 


248, 271 
373, 587 
265, 579 

2397 

2322, 2421 

1062, 1760 

2323, 2366 

1734, 2285 

1002 
996 
1008, 729 
716, 1001 
1.71 x 10-8 
1.39 x 1074 
0.002 


Table 1: All the masses in this table are in units of GeV. All the points satisfy the constraints 
described in Section 4. BMP1 and BMP2 are examples of Z-pole solutions, BMP3 and BMP4 
are representatives of Higgs-pole solutions, and BMP 5 is an example for Higgsino LSP solutions. 
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neutralino DM, while € = 1 for red and blue points. In this plot, the two dips around 45 
GeV and 62 GeV indicate the Z-pole and Higgs-pole solutions. Here, we want to make a 
comment that in focused scans, we also got points beyond the current LUX bounds but we 
have chopped them out and have displayed throughout this work only those solutions which 
are consistent with these bounds. By the way, if we introduce an axino as the LSP, i.e., the 
lightest neutralino is not the LSP, these chopped points are still natural and consistent with 
all the current experimental constraints. Moreover, we can see that in the near future the 
XENONIT experiment will completely probe solutions of our present scans. One can also 
notice that there is a wide gap between the Higgs-pole solutions and Higgsino LSP solutions 
(green points with mgo between 250 GeV 350 GeV). We notice that the o°!(x,p) is too high 
for points with neutralino mass between 65 GeV to 250 GeV. Even if we rescale the o% (x, p), 


points still rule out by the current LUX bounds. In addition to it, we also notice that for 


the Higgsino LSP scenario, x9, x}, and x7 are Higgsino dominated, x3 is Bino dominated, x? 


and yz are Wino dominated. Since mọ ~ myo ~ m,+, leptons from xg and y7 are hard 


x x 
to reconstruct. The most effective channels that could contribute to 3 leptons is pp > x{xy5. 


However, in this scenario, m o ~ X3 < 520 GeV, whereas the ATLAS Collaboration could only 
exclude points with Wino mass smaller than 380 GeV [56]. The CMS Collaboration has better 
results, but only excludes the points with Wino mass smaller than 450 GeV [57]. One can also 
see [58, 59] for probing light higgsino using monojet searches. This implies that in our case, 
even though the Higgsino LSP solutions are natural solutions but are out of reach of the present 
LHC searches. And we have confirmed it from numerical calculations for LHC SUSY searches 
as well. It is therefore, we will not consider them for further analyses. 

We want to comment on the light stau-neutralino coannihilation solutions. We find that 
if we insisting on Agw < 50, the light stau-neutralino coannihilation scenario is knocked out 
though it can be achieved if we relax Agw up to 100. This is why we will consider only Z-pole 
and Higgs-pole solutions for collider studies. 

We have collected five represented benchmark points (BMP) in Table 1. BMP1 and BMP2 
are the examples of Z-pole solutions with correct relic density and small relic density respec- 
tively. In these points, gluinos and the first two generations of squarks are heavier than 2 TeV, 
while top and bottom squarks tie and by have masses from 1 TeV to 2 TeV. The first two 
generations of sleptons have masses are in the range of [800, 1000] GeV while 72 are in [720, 
950] GeV mass range. For BMP1 and BMP2, Agw is 29 and 24 while Mgs is 333 GeV and 260 
GeV, respectively. BMP3 and BMP4 represent Higgs-pole solutions with correct relic density 


and small relic density, respectively. The colored sparticles have masses in same range as for 
BMP1 and BMP2. The first two generation sleptons are in the mass range [520, 790] GeV while 


10 
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72 are in [160, 790] GeV. BMP5 show a higgsino LSP solution. This point also has similar 
spectrum as BMP1 with relatively heavy sleptons and winos are heavier than 550 GeV. Since 
these points have very small relic density (Qh? ~ 0.002), we rescale the direct detection rate as 
£0°!(x,p). Moreover, we can see that BMP2, BMP3 and BMP4 have Aa, within 30 [36]. 


6 LHC searches 


In this Section, we examine the constraints from the direct searches for the SUSY particles at 
the LHC on samples with relic density consistent with or smaller than 54 WMAP9 bounds and 
also satisfy the current LUX limits on direct detection of LSP neutralino. In order to realize 
small fine-tuning and satisfy experimental observations simultaneously, only electroweakinos 
(neutralinos and charginos) and stau are light and could be explored at the current LHC. 


Therefore, we should consider the following electroweak Drell-Yan production processes: 


PPX, PPSA K pp UXT, (i,j = 2,3,4 and l,m = 1,2). (21) 
In general, the produced neutralinos could decay through on-shell or off-shell Z® or Ah): 
2. SE, (22) 


where the charginos could only decay through W“™), 


Xx >W y (23) 


When 7; is light, such as BMP 3 and BMP 4, new decay modes of ¥? and %7 are possible, 


Yom, and, XF > Fw, (24) 


and 7, decay into ¥? with branching ratio approximates to 100%. 
The relevance of different decay modes depends on mass spectrums and will significantly 


influence collider searches for these particles. In Fig. 2 we show branching ratios of ¥9 (Figs. 2(a) 


and 2(c)) and yf (Figs. 2(b) and 2(d)) for samples considered in this work. The dominant decay 
channel of 8 for samples of Z pole is ¥3 > ¥?Z“) when the mass difference msg — myo is small. 
Once the decay into Higgs boson is kinematically possible, branching ratio to Y?h increase with 
increasing of mo — myo and become the dominant channel when mo — myo Z 140 GeV. The 
decay channels of ¥f is always ¥f > ¥?W™). For samples of Higgs pole, situations are more 


complex due to light 7,, as can be seen from BMP3 and BMP4. The decay of ¥} to 7,7 would 


be significant or even dominant 2(c). 
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Figure 2: Decay branching ratios of ¥8 (a and b) and Xf (c and d) for samples with relic density 
consistent with or smaller than 54 WMAP9 bounds. 


ATLAS [60, 61] and CMS [57] have performed electroweakinos searches for wino ¥$/X7 
with particular decay models. We use the powerful package CheckMate [62, 63, 64] (where 
PYTHIA 8 [65] and a tuned version of DELPHES 3 [66] have been used internally) to implement 
LHC constraints. NLO production rates are obtained by rescaling LO rates with K-factors 
calculated by Prospino 2 [67], which yield about 1.2 for higgsino pair production. 

As for electroweakinos searches, currently CheckMate has only employed the ATLAS anal- 
yses with 13.3 fb’ data [60]. So in order to fully take into account the current constraints, we 
also recast the latest ATLAS [61] and CMS [57] analyses based on a Monte Carlo simulation. 
In the simulation, MadGraph 5 [68] is adopted to generate background and signal samples, and 
PYTHIA 6 [69] is employed to handle the parton shower, particle decay, and hadronization pro- 
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cesses. We use MLM scheme to deal with the matching between matrix element and parton 
shower calculations, and use Delphes 3 [66] to carry out a fast detection simulation with the 
CMS setup. Jets are reconstructed using the anti— kr algorithm [70] with a distance parameter 
AR = 0.4. 

Generally, heavy electroweakinos productions with successive decay will lead to multi- 
leptons signal, among them 31+ Ér and 2l + Ér give the best sensitivity at the LHC searches. 
In the case of 31 + Ér search channel, major SM backgrounds are ZZ and WZ productions. 
Two leptons from Z decay are required to form same-flavor-opposite-sign (SFOS) pair. Two 
useful kinematic variables to discriminate signals from backgrounds are mr and Ép, where 
mr is the transverse mass defined as mp = 4/2(p}p Ér — ph pms) with piss is the missing 


transverse momentum vector and the lepton / is the one not forming the SFOS lepton pair. 


pp collider, Vs = 13 TeV, 3 leptons pp collider, Vs = 13 TeV, 3 leptons 
10! gpg 10! pt 


events / bin fb! 


0 50 100 150 200 250 300 


E, (GeV) 
(a) mr distributions (b) Ér distributions 


Figure 3: mr (a) and Ær (b) distributions for backgrounds and signal benchmark points in the 
21+ Ér channel at the 13 TeV LHC. 


In Fig. 3 we present the mr and Ér distributions of backgrounds and signals. In the case of 
ZZ background, the 3l final state mainly comes from the decay of both Z boson into [*/~ pairs 
with one lepton do not be successfully reconstructed. As there is no neutrino contributing Ér, 
its Ér distribution is softer than others, and so is its mr distribution. For the WZ background, 
the mr variable is bounded by the W boson mass, leading to an obvious endpoint near mw. 
All mp and Éy for signals are harder and are easy to be distinguished from backgrounds. 

In the case of 21 + Éy search channel, dominant backgrounds are WZ, WW, ZZ, and tt 
production. Still, these two leptons from Z decay are required to form SFOS pairs, whose 


invariable is a useful variable to distinguish signals from SM backgrounds. Another useful 
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variable m7», is defined as 


mr = | min - {max[mr(pt, pr), Mr(Pr, PT)]}, (25) 
Ppt+Ppp=Pr 


where my(pi,, p>) = \/2(pi:p', — pi- p4), and p% and p? are the transverse momenta of two 


visible particles in the decay chain (two leptons in our case). pp and p% are a partition of 
the missing transverse momentum p's. By definition, Mrz is the minimum of the larger mr 
over all partitions, its distribution for two identical chains has an upper endpoint, which is 


determined by the mass difference between the parent particle and its invisible child. 


pp collider, Vs = 13 TeV, 2 leptons pp collider, Vs = 13 TeV, 2 leptons 
PERS EO LESS SGRS TREATS ORE RE E r 
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3 i =. i F kac ! --- tt —— BMP4 
= I 


events / bin fb” 


0 50 100 150 200 250 300 0 50 100 150 200 250 300 


my (GeV) M2 (GeV) 
(a) my distributions (b) mg distributions 


Figure 4: my (a) and mr, (b) distributions for backgrounds and signal benchmark points in 
the 2l + Ér channel at the 13 TeV LHC. 


In Fig. 4 we demonstrate the my, and m2 distributions of backgrounds and signals. For the 
WZ and ZZ backgrounds and signals, lepton pairs from Z boson decay result in peaks around 
mz in the my distributions, as shown in Fig. 4(a). Whereas these two leptons for WW and tt 
backgrounds origin from two particles and do not have obvious feature. Fortunately, the m2 
distributions for the WW and tt backgrounds are essentially bounds by my. 

For the analyses of Ref. [60], we use CheckMATE to calculate corresponding significance. And 
for the analyses of Refs. [57] and Ref. [61], we apply the same cuts in various signal regions, 
and compare the obtained cross sections to 95% limits tabulated in these literatures. In Fig. 5, 
we present the 95% C.L. exclusion results of the LHC electroweakino searches in the MoMo 
plane (Fig. 5(a)) and mgo-Arw plane (Fig. 5(b)). 31 + Ey searches require exactly three hard 
leptons. As a result, samples close to Z threshold, i.e., mgg ~ mgo + mz, are hard to explore. 
Points below this threshold would be excluded by 2l + Ér searches due to small Mt and large 


production cross sections. 
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Figure 5: 95% C.L. exclusion results of the LHC electroweakinos searches in the mgo-mgo plane 
(a) and myo-Apw plane (b). 


Roughly, 31+ Ér and 21+ Fp searches at the current LHC could exclude higgsino dominant 


X$/X7 with mass up to 300 GeV. This result is consistent with our previous prediction [71] 
whereas seems to somewhat weaker than that given by the ATLAS [60] and CMS Collabo- 


rations [57]. The main reason is that in their searches, pure wino ¥$/X7 is assumed, which 


have larger production cross sections. Besides, they assume %7 and $ decay via W* and Z* 
bosons with a branching fraction of 100%, whereas decay branching ratios in our samples highly 
depend on mass spectra. 

Another and even stricter constraint on samples of Higgs pole come from the searches for 


electroweakinos with tau final states. These searches have been performed by both CMS [57] 


and ATLAS [61] Collaborations, by latter it was shown that when ¥9/%7 decay into ¥? via an 


intermediate on-shell stau or tau sneutrino, ¥$/Y7 with mass up to 760 GeV are excluded for 
a massless ¥?. In our case samples with ¥$ mass up to 460 GeV for Higgs pole could still be 
excluded, as shown in Fig. 5(a). 

Finally, we project exclusion results into myo-Agw plane (Fig. 5(b)). Samples of Z pole 
with small Agw are easy to be explored, whereas these with large Agw are hard to be excluded 
due to large u, which in turn indicate large mo and small production cross sections. In the 
case of Higgs pole, many samples with Agw up to 50 could by excluded by electroweakinos 
searches with tau final states. For both Z pole and Higgs pole, samples with Agw approximate 


to 20 could still survive, indicating naturalness of this SUSY framework. 
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7 Discussions and Conclusion 


We have studied natural supersymmetry in the GmSUGRA, and found that after demanding 
Agw Á 50, only the parameter space related to Z-pole and Higgs-pole solutions are left. We 
performed the focused scans for such parameter space and showed that it satisfies various 
phenomenological constraints and is compatible with the current direct detection bound on 
neutralino DM reported by the LUX experiment. Such parameter space also has solutions with 
the correct DM relic density besides the solutions with relic density smaller or larger than 50 
WMAP9 bounds. We also performed the collider study of such solutions by implementing and 
comparing with relevant studies done by the ATLAS and CMS Collaborations. We showed 
that the points with the higgsino dominant ¥$/X;7 mass up to 300 GeV are excluded for Z pole 


scenario while for Higgs-pole scenario, points with ¥ mass up to 460 GeV are excluded. Next, 
we displayed that both for the Z-pole and Higgs-pole scenarios, the points having Agw ~ 20 
still survive. Moreover, we present five benchmark points as examples of our present scans. In 
these benchmark points, gluino and the first two generations of squarks are heavier than 2 TeV, 
top squarks fı, are in the mass range [1, 2] TeV, while sleptons are lighter than 1 TeV. We 
also discuss that stau-neutralino coannihilation scenario is not compatible with our demand of 
Agw Í 50. On the other hand higgsino LSP solutions which are natural solutions but are out 
of reach of present LHC searches. Some part of the parameter space can explain the anomaly 


of muon (g — 2), within 30 as well. 
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